Self-duality and bound states of the toric code model in a transverse field 
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We investigate tiie effect of a transverse magnetic field on the toric code model. We show that 
this problem can be mapped onto the Xu-Moore model and thus onto the quantum compass model 
which are known to be self-dual. We analyze the low-energy spectrum by means of perturbative 
continuous unitary transformations and determine accurately the energy gaps of various symmetry 
sectors. Our results are in very good agreement with exact diagonalization data for all values of the 
parameters except at the self-dual point where level crossings are responsible for a first order phase 
transition between a topological phase and a polarized phase. Interestingly, bound states of two and 
four quasiparticles with fermionic and bosonic statistics emerge, and display dispersion relations of 
reduced dimensionality. 
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Topologically ordered phases, such as those present 
in fractional quantum Hall systems^'"'*, have attracted 
much attention in the last few years. Indeed, in his sem- 
inal paper, Kitaev showed that topologically degenerate 
ground states may serve as a robust quantum memory \ 
while braiding of anyonic excitations can be used for 
fault-tolerant quantum computation. Topologically pro- 
tected qubits have left the realm of theory since super- 
conducting nanocircuits led to their first experimental re- 
alization' . Recent progress in the field of ultracold atoms 
trapped in optical lattices also promises an implementa- 
tion of such systems'"''"'^". 

Although a topological quantum memory is, by na- 
ture, protected from decohcrence, it is natural to wonder 
how large a local perturbation can be before this pro- 
tection fails. With respect to this problematics, the toric 
code model (TCM), which is undoubtedly one of the sim- 
plest model displaying topological order^ , is a perfect test 
ground. In the presence of parallel magnetic fields, the 
breakdown of the topological phase has been shown to be 
caused by single-anyon condcnsation^^d2,i3,i4^ leading to 
two second-order transition lines merging in a topological 
multicritical point '^ ''^ ^ 

The aim of this Rapid Communication is to investigate 
the infiuence of a transverse field in the TCM which turns 
out to display completely different physics. Indeed, as we 
shall see, this model can be mapped onto the self-dual 
Xu-Moore model proposed to describe superconducting 
arrays Note that the Xu-Moore model can also be 
mapped onto the quantum compass model^ ' relevant for 
orbitally frustrated systems and for topologically pro- 
tected qubits^'"*'^'^'^". All results given below are thus also 
valid for these two models as far as the energy spectrum is 
concerned. In the following, we compute the low-energy 
spectrum by means of perturbative continuous unitary 
transformations (PCUTs) and compare our results with 
exact diagonalization (ED) data. Our results reveal the 




FIG. 1: (Color online). Original square lattice on which 
plaquettes p and stars s are defined. Big red (small black) 
dots define the lattice A (A). Here, we show the lattice with 
N = 18 spins and (implicit) periodic boundary conditions. 
Contour Ci (C2) is one of the diagonal (anti-diagonal) cycles 
used to define conserved parity operators. 

existence of a first order phase transition at the self-dual 
point and emphasize the importance of strong binding 
effects leading to a plethora of multi-quasiparticle bound 
states with kinetics of reduced dimensionality. 

Model — The transverse-field TCM Hamiltonian reads 

as 

H = -jJ2As-jJ2B,-hyY,<^h (1) 

S p i 

where = ll^Gs'^i' = Uiep^f' ^nd the erf's are 
Pauli matrices. Subscripts s and p refer to stars (ver- 
tices) and plaquettes of a square lattice, whereas i runs 
over all bonds where spin degrees of freedom are located 
(see Fig. 1). In zero field, one recovers the TCM"' whose 
topological ground state has eigenvalue -1-1 for all As and 
Bp operators. Excitations are Z2-charges with eigenvalue 
— 1 for one As (Z2 -fluxes with eigenvalue —1 for one Bp) 
localized on the stars (plaquettes). These particles are 
hard-core bosons with mutual half-fermionic (semionic) 
statistics. Charges (or fluxes) can only appear in pairs for 
a system with periodic boundary conditions. In a mag- 
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netic field, elementary excitations become dressed any- 
onic quasi-particles (QP)'^ ^ In the opposite limit J = 0, 
the ground state is fully polarized. Elementary excita- 
tions are spin flips (magnons) , which are likewise dressed 
when switching on J. 

Although As and Bp are no longer conserved in a trans- 
verse field, the parity operator YiiiEC '^f ^^^^^ commutes 
with H, provided C is a diagonal or anti-diagonal contour 
such as the ones depicted in Fig. 1. In the s eigenbasis, 
the parities of the number of spin flips along such con- 
tours are thus conserved. This important property allows 
for ED of "rather large" systems up to iV = 32 spins, with 
the periodic boundary conditions defined in Ref. 4. The 
product of two parity operators defined on parallel con- 
tours is furthermore equal to the product of all Ag and 
Bp operators between the corresponding contours, which 
relates parities of magnons to that of anyons. 

Self-duality — This correspondence is only one signa- 
ture of the strong link between both types of QP which 
roots in a crucial property of the model : its self-duality. 
This feature directly stems from the mapping of the 
transverse-field TCM onto the Xu-Moore model which 
is self-dual'^ Indeed, let us introduce spin variables 

living on the dual lattice A (see Fig. 1) 

al=As, al=Bp, and SJ^H^"' (2) 

where js{p) denotes the center of a star (plaquette). 
The notation j > i defines the set of all sites i E A 
whose two coordinates are smaller than those of j G A. 
Hamiltonian (1) can then be rewritten as that of the Xu- 
Moore model 

where the first (second) sum is performed over all sites j 
(plaquettcs p) of A. Note that the above mapping only 
holds in the thermodynamic limit and for open boundary 
conditions, and that the infinite number of spins involved 
in the definition of ctJ cannot keep track of degeneracies. 
In particular, ED spectra of the transverse-field TCM 
with periodic boundary conditions discussed below are 
not symmetric under the exchange hy ^ J. 

Interestingly, Nussinov and Fradkin' ' showed that the 
Xu-Moore model can also be mapped onto the quan- 
tum compass model^'"". This model has focused much 
attention recently and latest numerical results plead in 
favor of a unique first order transition at the self-dual 
point'^^'^^'^'^ contrary to the original claim by Xu and 
Moore^'"''^'^'. This scenario that we shall confirm in the 
following immediately implies that the topological phase 
is rather well protected from transverse fields compared 
to parallel fields. Indeed, in the former case it breaks 
down at hy = J whereas in the latter case, the transition 
takes place for a field magnitude of order J/S^"*. 

Perturbative analysis — As in Ref. 14, a PCUT treat- 
ment can be set up in the limit of low (high) field, high- 
lighting the role of the corresponding QP. namely anyons 



-0.8 




FIG. 2; (Color online). Ground-state energy per spin eo ob- 
tained from PCUT and ED. The vertical dashed line marks 
the transition at 6* = tt/4. Inset: magnetization my — —dhyGo. 

(magnons). Basically, this method transforms H into 
an effective Hamiltonian unitarily equivalent to H but 
conserving the number of Qpi 1,24,25,26,27,28^ allows to 
investigate the thermodynamic limit and sorts the en- 
ergy levels according to their number of QP. We shall 
thus investigate low-energy sectors and confront the QP 
interpretation stemming from PCUT with ED spectra. 

OQP sector — By construction, the OQP state is the 
ground-state and lies in the symmetry subspace where 
all parities are even. We have computed the perturba- 
tive expansion of the ground-state energy per spin eo, up 
to order 10, in the low-field regime. Setting J = cos6', 
hy = sin 9, and t = tan 9, it reads as 

eo _ _^ t^ 13 1'^ 197 1^ 163885*8 
cos 61 ^ " " 8" ^ 1536 " 98304 " 226492416 
186734746441 t^" 
~ 1174136684544000' ^ ' 

This formula and its high-field counterpart (obtained by 
exchanging hy and J) are represented in Fig. 2. As can 
be seen, the agreement between (4) and ED results for 
= 32 spins is remarkable. Let us mention that the 
PCUT series expansion is rather well converged, since 
the difference between order 8 and 10 is of the order 10~* 
for all values of 9. Furthermore, a Pade approximant 
analysis gives eo{9 = 7r/4) = —0.8038(1), which perfectly 
matches previous numerical results'^"'-'. 

The cusp in the ground-state energy at 9 = tt/A in- 
dicates that the topological phase breaks down when 
hy reaches the value J, in agreement with the self- 
duality of the model. The transition point is best de- 
tected when looking at the magnetization in y-direction, 
which is obtained from the Hellmann-Feynman theorem : 
niy = —dhyeo- It displays a jump that reveals the first 
order nature of the transition (see inset of Fig. 2). Our 
perturbative treatment therefore confirms the order of 
the transition in the quantum compass model""'""' 

IQP sector — We now turn to the properties of a single 
QP. These excitations are static (dispersionless) due to 
parity conservation. In the high-field phase, they belong 
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FIG. 3: (Color online). Comparison between PCUT and ED 
results (A'' = 32) for the IQP gap Ai and lowest 2QP gaps 
Aj' ''^ ■ Pictograms give an illustration of the four correspond- 
ing states (with 7r/4-tilted lattice compared to Fig. 1). Crosses 
denote particles and filled circles empty sites. 



to sectors with all spin-flip parities even, except exactly 
one diagonal and one anti-diagonal parities, crossing at 
the QP's position. We computed the energy gap Ai of 
this IQP sector which, at order 10 and in the low-field 
regime, reads as 
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Once again, ED data perfectly match analytical results, 
as can be inferred from Fig. 3, whose upper left pictogram 
gives a representation of a IQP state. At the transition 
point hy = J, Pade extrapolations lead to Ai{9 — 7r/4) = 
0.9005(1). As already mentioned, eigenstates on both 
sides of the transition have to be interpreted either in 
terms of dressed anyons or in terms of dressed magnons. 
Furthermore, we remark that ED, performed on clusters 
with periodic boundary conditions, can only detect the 
IQP excitation at large field since the excitation of a 
single anyon in the topological phase is forbidden for such 
boundary conditions. We also note that the one-magnon 
state is connected to one of the two-anyon states which 
we study below. 

2QP sector — In the high-field phase, parity symme- 
tries of a two-magnon state can be of two kinds. The 
first is obtained by setting all parities even except two 
diagonal or anti-diagonal parities which are odd as in 
the two lowest pictograms in Fig. 3. In this case, the 
two magnons can only move in the direction orthogo- 
nal to their relative position. This is a nice illustration 
of the dimensional reduction phenomenon ^"''i*' in which 
the transverse magnetic field induces a one-dimensional 
correlated hopping and leads to the formation of bound 
states. The second kind is obtained by setting all pari- 
ties even except two diagonal and two anti-diagonal par- 
ities which are odd, as depicted for instance in the upper 
right pictogram in Fig. 3. In such a configuration, par- 



ity conservation imposes very limited kinetics of the two 
magnons. Again, the transverse magnetic field leads to 
strong binding effects. 

The self-duality allows for a similar analysis in the 
topological phase for anyons, with the restriction that 
only two-charge or two-fiux states are allowed for a sys- 
tem with periodic boundary conditions. Such excita- 
tions have bosonic statistics. However, dyonic bound 
states made of one charge and one flux (with fermionic 
statistics) only exist for open boundary conditions. We 
have calculated all 2QP excitation energies up to order 
8. Hereafter, we provide the three lowest gaps in the low- 
field phase corresponding to the three 2QP configurations 
shown in the pictograms of Fig. 3, 
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These PCUT scries, plotted in Fig. 3, seem to indi- 
cate that Ai and A| are equal at the transition. As 
previously, we compared them with ED data which, for 
A'^ = 32 spins, imply to deal with blocks containing up to 
16 million states. Although ED and PCUT results almost 
lie on top of each other, ED spectra reveal the formation 
of energy jumps at the transition point, which can only 
be explained by level crossings with higher-energy states 
occurring in the thermodynamical limit. These cross- 
ings cannot be captured by the PCUT approach, whose 
perturbative nature imposes an adiabatic continuation of 
levels. However, we insist on the validity of our results 
for all 0's, except at the transition point. 

We therefore conclude that the level crossing respon- 
sible for the cusp in the ground-state energy does not 
originate from IQP and 2QP levels, since these excita- 
tion energies are finite at 9 — it/ A. 

4QP sector — To address the origin of the cusp, we 
now look for the lowest excited state belonging to the 
same symmetry sector as the ground state. PCUT en- 
ergy ordering suggests a 4QP state as a natural candi- 
date. Such a 4-magnon state (or a two-flux and two- 
charge state with bosonic statistics) is built from all con- 
figurations where the magnons occupy the corners of a 
rectangle. These can be linked to the configuration where 
the four QP form a close-packed square, by shifting the 
center of mass and/or the relative positions of the QP. 
Four such configurations are shown in Fig. 4. 

In contrast to 2QP states, 4QP states in this parity sec- 
tor have a two-dimensional dispersion. However, a partial 
dimensional reduction still occurs for the relative motion 
of the QP. Indeed, the corresponding effective Hamilto- 



4 




PCUT 
order 4 



FIG. 4: (Color online). Lowest 4QP gap A4 obtained from 
PCUT and ED. Pictograms illustrate four 4QP configura- 
tions, with different relative positions. 



nian at order n in perturbation is found to be that of 
a single particle moving in n coupled one-dimensional 
chains, with an impurity whose extension grows with n 
(details will be given elsewhere). As it frequently occurs 
in this type of problem, the bound state associated to this 
imurity cannot be obtained perturbatively, and one has 
to resort to a numerical diagonalization of the effective 
Hamiltonian. 

The gap A4 of the 4QP bound state, obtained from the 



fourth-order effective Hamiltonian, is shown in Fig. 4, to- 
gether with ED results. Both match away from the tran- 
sition but, contrary to the ED gap which goes to zero 
at the transition, the PCUT gap remains finite. At or- 
der 3, one gets A4 = 1.728 at 6' = it/ A whereas at order 
4, A4 = 1.721 suggesting a fast convergence. This fi- 
nite value shows once again that PCUT miss level cross- 
ings occurring at the self-dual point, but give reliable 
results everywhere else. This discrepancy can be readily 
explained by the first order transition, in which a cas- 
cade of level crossings from high-energy states down to 
the ground state occurs in the thermodynamic limit. 

Perspectives — We have shown that a transverse mag- 
netic field in the TCM is the source of important bind- 
ing effects, leading to a sequence of bound states with 
reduced dimensional kinetics, in deep contrast with the 
parallel field case' '. The fate of these bound states in 
a general magnetic field, where single-quasiparticle exci- 
tations are also dispersive, is a fascinating issue left for 
future studies. 
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